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Abstract The potential energy surfaces of even-even 142_156 Ba are investigated in the constrained reflection- 
asymmetric relativistic mean-field approach with parameter set PKT It is shown that for the ground states, 
142 Ba is near spherical, 156 Ba well quadrupole-deformed, and in between 144_154 Ba octupole deformed. In 
particular, the nuclei 148,150 Ba with N = 92,94 have the largest octupole deformations. By including the 
octupole degree of freedom, energy gaps N — 88, TV = 94 and Z = 56 near Fermi surfaces for the single- 
particle levels in 148 Ba with fa ~0.26 and fa ~0.17 are found. Furthermore, the performance of the octupole 
deformation driving pairs (^2/ 7/ / 2 , vli 13 / 2 ) and (n2d 5 / 2 , is demonstrated by analyzing the single- 

particle levels near Fermi surfaces in 148 Ba. 
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' 1 Introduction 

00 

00 

In the last decades, the phenomena related to octupole deformation have received wide attention. Normally 
the regions of nuclei with strong octupole correlations correspond to either the proton or neutron numbers close 
OO 1 to 34 (lg 9/ 2 <-> 2p 3/2 coupling), 56 (lh 11/2 O 2d 5/2 coupling), 88 (li 13/2 2/ 7/2 coupling), and 134 (lj 15 / 2 2sg /2 
^ , coupling) [l|. 

Extensive efforts have been made for understanding the structure of neutron-rich nuclei around Z ~ 56 and 
J* ■ N ~ 88. Experimentally in this region, many octupole deformed bands have been identified and extended to 
• Ch ! higher spin, such as in 139 Xe [IS], in even-even 140 " 148 Ba in 144 < 146 Ce Q, and in 145 - 147 La Q. On the 

' theoretical side, a variety of approaches have been applied to investigate the role of octupole degree of freedom 
■ in this nuclear region. The Woods-Saxon-Bogoliubov cranking model is used to study the shapes of rotating 
Xe, Ba, Ce, Nd, and Sm nuclei with N = 84 — 94 and the expectations of octupole-deformed mean- fields at low 



and medium spins are confirmed lOj. Based on the Hartree-Fock plus BCS approximation and adiabatic time- 



dependent Hartree-Fock plus zero point energy in the cranking approximation, the energy splitting and B{E1) 



transition are well described for 140 Ba, 142 ~ 150 Ce, 144 " 152 Nd, and 14fi - 154 Sm [ll|. In the spdf interacting boson 
model, good agreement of the calculated low-lying energy spectra and transition rates with data is obtained 
for 144 Ba and 146 Ba [l3|. The reflection-asymmetric shell model is applied to describe the octupole deformed 
bands in neutron-rich 142 Ba and 145 Ba, and good agreement with the ex peri mental data is obtained (lj| . 



During the past years, the Relativistic Mean-Field (RMF) theory 14H16] has achieved great success in 



describing many nuclear phenomena related to stable nuclei 14j , exotic nuclei 17J, Il8| as well as supernova 
and neutron stars [3]. Recently, the Reflection- ASymmetric Relativistic Mean-Field (RAS-RMF) approach 
considering the octupole degree of freedom is developed and applied to the well-known octupole deformed 



nucleus 226 Ra 20], and La isotopes 2l|. In Ref. [22j, the RAS-RMF approach has been applied to investigate 
the Potential Energy Surfaces (PESs) of even-even 146_156 Sm isotopes in the (/3 2 , fiz) plane, and it is suggested 
that the critical-point candidate nucleus 152 Sm marks the shape/phase transition not only from U(5) to SU(3) 
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symmetry, but also from the octupole deformed ground state in 150 Sm to quadrupole deformed ground state in 
154 Sm. Therefore, it is interesting to investigate the Ba isotopes in the RAS-RMF approach. 

In this paper, the RAS-RMF approach will be applied to investigate the potential energy surfaces of even- 
even 142_156 Ba isotopes in the (/3 2 , Pz) plane, and the single-particle levels near Fermi surfaces for the nucleus 
148 Ba will be studied. 

2 Formalism 



The basic ansatz of the RMF theory is a Lagrangian density where nucleons are described as Dirac particles 
which interact via the exchange of various mesons and the photon. The mesons considered are the isoscalar- 
scalar a, the isoscalar- vector uj and the isovector- vector p. The effective Lagrangian density reads 23[, 



C 



i^dfi -M — g a a - g^u^ 

1-T 3 



~9p1 T ' Pfi e 7 



(1) 



in which the held tensors for the vector mesons and the photon are respectively defined as, 



(2) 



F„ 



Using the classical variational principle, one can obtain the Dirac equation for the nucleons and the Klein- 
Gordon equations for the mesons. To solve these equations, we employ the basis expansion method, which 
has been widely used in both the non-relativistic and relativistic mean-held models. For axial-symmetric 
reflection-asymmetric systems, the spinors are expanded in terms of the eigenfunctions of the Two-Center 
Harmonic-Oscillator (TCHO) potential 



f 



V(r ± ,z) = -Mu;iri 



(3) 



-MujI{z~z 2 )\ 



z^O 



where M is the nucleon mass, z 1 and z 2 (real, positive) represent the distances between the centers of the 
spheroids and their intersection plane, and Wi(o;2) are the corresponding oscillator frequencies for z < (z ^ 
0) [2^]. The TCHO basis has been widely used in the studies of fission, fusion, heavy- ion emission, and various 
cluster phenomena 24| . By setting proper asymmetric parameters, the major and the z-axis quantum numbers 
are real numbers very close to integers, and the integers are used in the Nilsson-like notation Q[Nn z mi] for 
convenience. More details can be found in Ref. 20| . 

The binding energy at a certain deformation is obtained by constraining the mass quadrupole moment (Q 2 ) 
to a given value fj, 2 [25[, i.e., 

(H') = {H) + \c({Q 2 )-n 2 f (4) 

where C is the curvature constant parameter, and p. 2 is the given quadrupole moment. The expectation value 
of Q 2 is (Q 2 ) = {Q 2 ) n + {Q 2 ) P with (Q 2 )„ yP = (2r 2 P 2 (cos#))„ iP . The deformation parameter /3 2 is related to (Q 2 ) 
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by, (Q 2 ) = Ar 2 /3 2 with r = i?, A 1//3 (R — 1.2 fm) and A the mass number. The octupole moment constraint 

V57T 

3 

can also be applied similarly with (Q 3 ) = (Q 3 ) n + (Q 3 ) P , (Q 3 )n, P = (2r 3 P 3 (cos6»)) n>p , and (Q 3 ) = —=Ar 3 j3 3 . By 

V77T 

constraining the quadruple moment and octupole moment simultaneously, the total energy surface in (/3 2 , $3) 
plane can be obtained. 

3 Results and discussion 

The properties of even-even 142_156 Ba are calculated in the constrained RAS-RMF approach with parameter 
set PK1 [26l | . The TCHO basis with 16 major shells for both fermions and bosons is used. The pairing 
correlation is treated by the BCS approximation with a constant pairing gap A = 11.2/yA MeV. 

To investigate the shape evolution in 142_156 Ba, the total energies as functions of 02 and f3 3 have been 
analyzed. As an example, Fig. Q] displays the contour plots for 142 Ba and 148 Ba. It is found that for the 
ground states, 142 Ba is near spherical without octupole deformation, 144_154 Ba octupole deformed and 15B Ba 
well quadrupole-deformed. In detail, for 144_156 Ba, the quadrupole deformation (3 2 of the global octupole 
minimum gradually increases with increasing neutron numbers. On the other hand, the octupole deformation 
j3 3 of the global octupole minimum gradually increases for 144,146 ' 148 Ba, and decreases for 15 °. 152 > 154 Ba, and until 
156 Ba the global minimum is well deformed with little octupole deformation. A soft area covering the global 
octupole minimum and the saddle point at /3 3 =0 appears in 144 Ba, and develops in 146_154 Ba. Furthermore, 
the energy difference between the global octupole minimum and the saddle point increases from 144 Ba to 148 Ba, 
while it decreases from 150 Ba to 156 Ba. 




Fig. 1. (Color online) The contour plots of total energies for even-even 142 Ba and 148 Ba in (/3 2 , f3 3 ) plane 
obtained in RAS-RMF approach with PK1 and constant-A pairing. The energy separation between contour 
lines is 0.5 MeV. The global minimum and other local minima are denoted by "•" and "▼" respectively. 

The binding energy, quadrupole and octupole deformations are summarized for the ground states of even- 
even 142_156 Ba in TableQ] The binding energies are well reproduced within 0.3%. Moreover, excellent agreement 
is obtained for the quadrupole deformations. It's indicated that the ground states of even-even 144_154 Ba with 
N = 88— 98 are octupole deformed while in the middle the nuclei 148 ' 150 Ba with N = 92 — 94 are the most 
octupole deformed. It is noted that the suggested octupole Ba nuclei (N = 88 — 98) are more neutron-rich than 
those (^ = 86 — 91) predicted in the finite-range droplet model [2^]. This conclusion consists with the previous 
RAS-RMF calculation for La isotopes (ilj . 

Table 1. The total binding energy (in MeV) as well as the quadrupole deformation /3 2 and octupole deforma- 
tion fi 3 of the ground states of even-even 142_156 Ba obtained in the constrained RAS-RMF approach with 
PK1, in comparison with the available experimental data. 



No. X 



ZHANG Wei et al: Octupole deformation for Ba isotopes in a 
reflection-asymmetric relativistic mean-field approach 



■1 



Nucleus 


Jjjca.1 


/3= al 


,3§ al 


gexp [27] 


/3^ xp [28] 


1 AO T-i 


1181.25 


0.12 


0.00 


1180.14 


0.16 


144 Ba 


1190.19 


0.20 


0.12 


1190.23 


0.19 


146 Ba 


1199.06 


0.23 


0.13 


1199.60 


0.22 


148 Ba 


1207.26 


0.26 


0.17 


1208.76 




150 Ba 


1215.12 


0.28 


0.17 


1217.55 




152 Ba 


1221.86 


0.30 


0.12 


1225.58 




154 Ba 


1228.30 


0.32 


0.08 






156 Ba 


1234.58 


0.33 


0.03 







To understand the evolution of the octupole deformation, the microscopic single-particle levels are analyzed. 
Fig. [5] shows the neutron single-particle levels for the states minimized with respect to f3 3 and the states with 
/3 3 =0 for (3 2 = 0.14 ~ 0.38 in 148 Ba side by side. The levels near the Fermi surface are labeled by Nilsson- 
like notations fl[Nn z mi] of the largest component. In the left panel of Fig. [2j an energy gap TV = 88 fenced 
with several levels near the Fermi surface for f3 2 = 0.20 ~ 0.30 is found approaching the octupole minimum. 
Additionally, a gap at N = 94 is found for j3 2 ~ 0.28. For the states with ^3=0, the neutron gap at N = 88 
is relatively small, and the gap at N = 94 does not appear. Furthermore, the proton single-particle levels for 
the corresponding states are shown in Fig. [3l An energy gap Z = 56 appears for the states with octupole 
deformation approaching the ground state in the left panel of Fig. [3l while the energy gap Z — 56 can not 
be found near the Fermi surfaces for the states with [3 3 =0 in the right panel of Fig. [3] Thus the energy gaps 
N = 88, N = 94, and Z — 56 near Fermi surfaces are responsible for the octupole minimum of Ba isotopes. 
In particular, the energy gaps TV = 94 and Z = 56 around (3 2 ~ 0.28 are consistent with the large octupole 
deformation predicted for 148 ' 150 Ba. 




0.16 0.20 0.24 0.28 0.32 0.36 0.16 0.20 0.24 0.28 0.32 0.36 



4 4 

Fig. 2. (Color online) Neutron single-particle levels of 148 Ba in RAS-RMF approach with PK1 as functions 
of @2 for states minimized with respect to fiz (left panel) and states with /?3 = (right panel). The dash-dot 
lines denote the corresponding Fermi surfaces. The levels near the Fermi surface are labeled by Nilsson-likc 
notations Q,[Nn z mi] of the first component. The corresponding /J3 are shown in the inset. The quadrupole 
deformation of the ground state is indicated by the vertical gray line. 
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0.16 0.20 0.24 0.28 0.32 0.36 0.16 0.20 0.24 0.28 0.32 0.36 

Fig. 3. (Color online) Same as Fig. but for proton. 

It is well-known that for nuclei with N ~ 88 or Z ~ 56 the octupole deformation driving pairs of orbitals 
include (i/2/ 7 / 2 , ^1*13/2) and (7r2d 5/ / 2 , 7rl/iii/ 2 ), which in the axially deformed case will be subgrouped as 
(i/l/2[541], 1/1/2 [660]), (i/3/2[532], z/3/2[651]), (i/5/2[523],i/5/2[642]), (i/7/2[514], 1/7/2 [633]), and (tt1/2[431], 
7rl/2[550]), (tt3/2[422], tt3/2[541]), (tt5/2[413], tt5/2[532]), respectively. It is interesting to investigate the 
performance of such pairs in the single-particle levels near Fermi surfaces in Figs. [2] and El These levels 
together with their BCS occupation probabilities and their ten leading components for the ground state are 
shown in Tables [5] and [3] respectively. Taking the level i/3/2[521] as an example, its second (18.7%) and 
third (16.5%) components compose an octupole deformation driving pair (i/3/2[532], i/3/2[651]). Similarly, 
one can find the pair (1/2[541], 1/2 [660]) for z/l/2[530], the pair (3/2[532], 3/2[651]) for i/3/2[532], and the 
pair (5/2 [523], 5/2 [642]) for z/5/2[523]. For the proton side, octupole deformation driving pairs are also found 
among the components of single-particle levels near Fermi surfaces at the ground state of 148 Ba. Therefore 
both the neutron and the proton driving pairs play important roles for the octupole minimum in 148 Ba. 

Table 2. Single-particle levels near Fermi surface for the ground state in 148 Ba together with their BCS 
occupation probabilities and corresponding contributions from the ten leading components. The components 
originating from the octupole deformation driving pairs of orbitals (u2fj/ 2 , VU13/2) and (7r2d 5 / 2 , Tl/111/2) 
are in bold. 



level 


yl/2[530 




i/3/2[532 




iv3/2[521] 


1/5/2 [523 




occu. 


0.973 




0.968 




0.721 




0.609 




1st comp. 


l/2[530] 


41.0% 


3/2[532] 


52.6% 


3/2[521] 


36.0% 


5/2[523] 


56.4% 


2nd comp. 


1/2 [541] 


15.6% 


3/2[541] 


17.6% 


3/2[532] 


18.7% 


5/2[532] 


12.7% 


3rd comp. 


l/2[660] 


6.2% 


3/2[512] 


6.2% 


3/2[651] 


16.5% 


5/2[642] 


10.1% 


4th comp. 


1/2[651] 


5.8% 


3/2[651] 


6.0% 


3/2[631] 


6.8% 


5/2 [633] 


5.0% 


5th comp. 


1/2[510] 


5.7% 


3/2[631] 


3.7% 


3/2[642] 


6.7% 


5/2[622] 


3.3% 


6th comp. 


1/2[631] 


3.5% 


3/2[521] 


2.2% 


3/2[761] 


2.5% 


5/2 [503] 


3.2% 


7th comp. 


l/2[640] 


3.2% 


3/2[761] 


1.7% 


3/2[501] 


2.1% 


5/2[512] 


2.0% 


8th comp. 


l/2[770] 


3.1% 


3/2[402] 


1.6% 


3/2[622] 


2.0% 


5/2 [413] 


1.3% 


9th comp. 


l/2[400] 


2.6% 


3/2[642] 


1.1% 


3/2[721] 


0.9% 


5/2[752] 


0.8% 


10th comp. 


1/2[521] 


1.9% 


3/2[321] 


0.9% 


3/2[752] 


0.8% 


5/2[613] 


0.5% 



Table 3. Same as Table [21 but for proton. 
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level 


?r3/2[422 




tt1/2[420; 




7r5/2[4l3 




7t3/2[541 




occu. 


0.936 




0.863 




0.098 




0.068 




1st comp. 


3/2 [422] 


69.4% 


l/2[420] 


38.6% 


5/2[413] 


78.7% 


3/2[541] 


40.8% 


2nd comp. 


3/2[431] 


10.7% 


l/2[550] 


14.4% 


5/2 [422] 


8.9% 


3/2[521] 


15.3% 


3rd comp. 


3/2 [402] 


5.1% 


l/2[530] 


11.0% 


5/2 [523] 


4.2% 


3/2[411] 


13.9% 


4th comp. 


3/2[541] 


3.8% 


1/2[541] 


8.6% 


5/2[532] 


2.0% 


3/2 [532] 


10.2% 


5th comp. 


3/2 [532] 


2.0% 


l/2[400] 


4.2% 


5/2 [303] 


1.7% 


3/2[422] 


7.1% 


6th comp. 


3/2[521] 


1.7% 


l/2[521] 


3.7% 


5/2[512] 


0.5% 


3/2[651] 


2.2% 


7th comp. 


3/2[312] 


1.4% 


l/2[660] 


3.3% 


5/2[813] 


0.1% 


3/2[512] 


1.9% 


8th comp. 


3/2[211] 


0.7% 


1/2[431] 


2.7% 


5/2 [633] 


0.1% 


3/2 [431] 


1.1% 


9th comp. 


3/2[651] 


0.5% 


l/2[220] 


1.6% 


5/2 [833] 


0.1% 


3/2[761] 


0.9% 


10th comp. 


3/2 [642] 


0.3% 


1/2[411] 


1.5% 


5/2[312] 


0.1% 


3/2[321] 


0.6% 



4 Summary 

In conclusion, the PESs of even-even 142_156 Ba in ((3 2 , (3 3 ) plane are investigated by the constrained RAS- 
RMF approach, and the single-particle levels near Fermi surfaces for the nucleus 148 Ba are studied. It is shown 
that for the ground states, 142 Ba is near spherical without octupole deformation, 144 ~ 154 Ba octupole deformed 
and 156 Ba well quadrupole-deformed. The nuclei with largest octupole deformation /? 3 for the ground states 
are predicted to be 148 > 150 Ba (N — 92,94). 

By including the octupole degree of freedom, energy gaps N — 88, N = 94, and Z — 56 near Fermi surfaces 
are found in single-particle levels approaching the ground state of 148 Ba. Furthermore, the performance of the 
octupole deformation driving pairs (v2f 7 / 2 , v^in/2) and {'K2d s /2 1 is demonstrated by analyzing the 

components of the single-particle levels near Fermi surfaces. 

The authors gratefully acknowledge Professor MENG Jie, Dr. GENG Li-Sheng and LIANG Rao-Zhao for 
their helpful suggestions and stimulating discussions. 
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